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ABSTRACT 

We consider ion heating by turbulent Alfven waves (AWs) and kinetic Alfven waves (KAWs) with wave- 
lengths (measured perpendicular to the magnetic field) that are comparable to the ion gyroradius and frequen- 
cies CO smaller than the ion cyclotron frequency Q.. As in previous studies, we find that when the turbulence 
amplitude exceeds a certain threshold, an ion's orbit becomes chaotic. The ion then interacts stochastically with 
the time-varying electrostatic potential, and the ion's energy undergoes a random walk. Using phenomenolog- 
ical arguments, we derive an analytic expression for the rates at which different ion species are heated, which 
we test by simulating test particles interacting with a spectrum of randomly phased AWs and KAWs. We find 
that the stochastic heating rate depends sensitively on the quantity e ~ 5vp /v±, where (vjj) is the component 
of the ion velocity perpendicular (parallel) to the background magnetic field Bq, and 5vp (5Z?p) is the rms ampli- 
tude of the velocity (magnetic-field) fluctuations at the gyroradius scale. In the case of thermal protons, when 
e <C Eciit, where Ecrit is a dimensionless constant, a proton's magnetic moment is nearly conserved and stochas- 
tic heating is extremely weak. However, when e > Ecrit, the proton heating rate exceeds the cascade power that 
would be present in strong balanced KAW turbulence with the same value of 5vp, and magnetic-moment con- 
servation is violated even when co <C il. For the random-phase waves in our test-particle simulations, Ecrit — 0.2. 
For protons in low-p plasmas, E 2± p^^/^5Bp/Z?0i and E can exceed Ecrit even when 5Bp/Bo ^ Gcrit^ where p is 
the ratio of plasma pressure to magnetic pressure. The heating is anisotropic, increasing much more than v| 

when p <C L (In contrast, at P > 1 Landau damping and transit-time damping of KAWs lead to strong parallel 
heating of protons.) At comparable temperatures, alpha particles and minor ions have larger values of E than 
protons and are heated more efficiently as a result. We discuss the impUcations of our results for ion heating in 
coronal holes and the solar wind. 

Subject headings: solar wind — Sun: corona — turbulence — waves — MHD 



L INTRODUCTION 

Beginning in the 1960s, a number of authors developed 
steady-state hydrodynamic models of the solar wind, in which 
the temperature was fixed at the coronal base and the solar 
wind was heated by thermal conduction (e.g. Parker 1965; 
Hartle & StuiTock 1968; Durney 1972; Holzer & Leer 1980). 
For realistic values of the coronal temperature and density, 
these models were unable to reproduce the large flow ve- 
locities of fast-solar-wind streams at 1 AU, suggesting that 
the fast wind is heated above the coronal base by some ad- 
ditional mechanism. Observational evidence for extended, 
non-conductive heating has since been provided by mea- 
surements from the Ultraviolet Coronagraph Spectrometer 
(UVCS), which show radially increasing minor-ion temper- 
atures in coronal holes (the open-magnetic-field-line regions 
from which the fast wind emanates) at heliocentric distances r 
between I.SRq and 3.5Rq (Kohl et al. 1998; Antonucci et al. 
[2000). Identifying the physical mechanisms responsible for 
this heating and determining the heating rates of the different 
particle species are among the major challenges in the study 
of the solar wind at the present time. 

One of the first mechanisms proposed to account for 
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solar-wind heating was turbulence (Coleman) 119681) . The 
importance of turbulent heating is suggested by in situ 
measurements of ubiquitous, large-amplitude fluctuations 
in the velocity and ma gne ti c fiel d in the interpl a netary 
medium ([B elcher and Davii 119711; iGoldstein elaLl 1 19951 
iBruno andC arbone 2005), as well as the positive correlation 
between the solar-wind temperature and the amplitude of the 
fluctuations dGrappin et al.il 1990HVasquez et al.ll2007l) . In ad- 
dition, the expected rate at which the measured fluctuations 
dissipate (based on phenomenological turbulence theories) is 
compara ble to th e observation ally infer red solar-win d heat- 
ing ra te (ISmith e t al. 2001; Bre ech et all 2009; Cranme retaTI 
120091) . The in situ measurements on which the above studies 
are based are limited to the locations where spacecraft have 
flown - that is, to r > 0.3 AU. However, the velocity and 
magnetic-field fluctuations are often correlated in the sense 
of Alfven waves propagati ng away from the Sun in the solar- 
wind frame (Belcher and DavisI [l97Tt iTu and Ma rs"chlll995l; 
^avassano et al. 2000), indicating that these waves originate 
at or near the Sun, consistent with the idea that turbulent heat- 
ing remains important as r decreases below 0.3 AU. 

At least two different scenarios for turbulent heating of 
coronal holes and the solar wind are possible. In the 
first, magnetic reconnection or some other process launches 
Alfven waves into the corona, including waves with > 
k±, where k^^ and k± are the components of the wavevec- 
tor k parallel and perpendicular to the local background 
magnetic field Bq.^ Given the large Alfven speed va in 

* Alfven waves play a key role in extended-heating models because fast 
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coronal holes (> 10^ km/s), the frequency co — fe||VA of 
such waves exceeds 1 Hz for wavelengths shorter than ^ 
10^ km. Once such waves enter the corona, nonlinear in- 
teractions with coronal density fluctuations [which are in - 
ferred from radio observations (IColes and HarmonI Il989l) 1 
can convert a significant fractio n of the A l fven w ave power 
into fast magnetosonic waves (iChandranI 120081) . The en- 
ergy in fast magnetosonic waves can then cascade to higher 
frequencies (Cho and Lazarian 2002; Svidzinski et al. 201Cj), 
generating high-frequency Alfven waves with > k± in 
the process (IChandranI l2005h . Although the dissipation of 
high-frequency fast waves and Alfven/ion-cyclotron waves 
could in principle explain the UVCS observations of ion 
heating ([Li and Habbal 2001; Hollweg and Isenberg 2002; 
iMarkovskii et al.ll2010l) . there is no direct observational evi- 
dence that waves with high frequencies and/or | > k± are 
present in coronal holes. 

An alternative scenario, which we focus on in this pa- 
per, involves the launching of much lower-frequency Alfven 
waves by convective photospheric motions. An effective 
k± for such waves can be estimated as 2n/Lo, where Lq 
is of order the average spacing of either supergranules 
(~ 3 X 10"* km) or photospheric flux tubes 5000 km). 
For a wave period 27t/(fc||VA) of 10^ s and an Alfven 
speed of 10^ km/s, the ratio \k^^\/k± of such waves is < 
0.03. Such highly anisotropic Alfven waves are ineffi- 
cient a^^generating compressive modes ( Cho and Lazaria^ 
120031; IChandranI 120051. |2008). On the other hand, they can 
interact with oppositely propagating Alfven waves, caus- 
ing wave energy to cascade from large scales to small 
scales, or, equivalently, small k to large k, whe re the fluc- 
tuations dissipate, heat ing the ambient plasma (llroshnikovl 
119631: iKraic hnan 1965). Although the Sun launches only 
outward-propagating Alfven waves, the inward-propagating 
waves required for the Alfven-wave cascade are generated 
near the Sun by non-WKLB wave reflection arising from 
the gradient in the Alfven spee d ( Heinemann and Olbe^ 
1980; 'Vein et al."1989'; 'Matthaeus et al. 1999; DmitruketaT* 
2002; Cranmer and van Ballegooijen 2005; Verdini and Velli 
20071: iHollweg and Isenbergll2007l: IVerdini et al.l l2009fl). An 
important development in the theory of Alfven-wave tur- 
bulence was the discovery that interactions between op- 
positely propagating Alfven waves cause wave energy 
to cascade primarily to larger k \ and only weakly to 
larger \k\\\ ('Montgomery and Turner 1981; Shebali n et al.l 
[1983:; .Goldreich and Sri dhar 1995 ). At the dissipation scale, 
the value of |^|||/^x is thus even smaller than at the driving 
scale Lq. 

This second scenario for solar-wind heating is compelling 
for a number of reasons. For example, convective pho- 
tospheric motions inevitably launch low-frequency Alfven 
waves into the corona by perturbing the footpoints of open 
magnetic field lines, and low-fre quency Alfve n waves are 
observe d in the corona ( Tomczvk et al.l 120071) and at r > 
0.3 AU ("B elcher and Davislll971i) . In addition, several mod- 
els have been developed to describe wave reflection and 
turbulent heating by low-frequency Alfven waves in the 

magnetosonic waves can not in general escape from tlie chromosphere into 
the corona since they are reflected at the transition region I Hollweg 1978). In 
addition, slow magnetosonic waves are strongly damped in coUisionless low- 
p plasmas ( Barnes 1966) and thus are not an effective vehicle for transporting 
energy from the coronal base to r > l.5R@. 



fast solar wind, taking into account the solar-wind veloc- 
ity, density, and magnetic-field profiles, and incorporating 
observational constraints on the Alfven-wave amplitudes; in 
all of these models, the turbulent heating rate appears to 
be consistent with the requirements for generating t he fast 
wind (Cranmer and van Ballegooijen 2005; Cranmer etal J 
2007; Chandran and Hollweg" 2009; Verdi ni et all I2009M 
i2010i) . We also note that radio observations of density fluctu- 
ations provide an upper limit on the heating rate from Alfven 
waves in coronal holes, since the Alfven waves become in- 
creasingly compressive with increasing k. Although these up- 
per limits rule out fast-wind generation by non-turbulent high- 
frequency Alfven/ion-cyclotron waves (unless k is nearly par- 
allel to Bo for all the waves; Hollweg 2000), they are con- 
sistent with fast-wind generation by low- frequenc y (kinetic) 
Alfven-wave turbulence with k± ^ k^^ (iHarmon and Coleg 
2005; Chandran et al. 2009). 

Despite these considerations, it is not clear that low- 
frequency Alfven-wave turbulence can explain two key obser- 
vations. First, measurements of the proton and electron tem- 
perature profiles in the fast solar wind at r > 0.3 AU demon- 
strate that the proton heating rate exceeds the electron heating 
rate by a modest factor (Cranmer et al. 2009). Similarly, em- 
pirically constrained fluid models of coronal holes including 
thermal conduction suggest that protons receive a substan- 
tial fr action (^ 0.5) of the total heating power (lAUen et al.l 
Il998h . Second, UVCS observations show that minor ions 
such as 0+^ are heated in such a way that thermal mo- 
tions perpendicular to Bq are much mor e rapid than ther- 
mal motions along B p (i.e., 71 > 7||) (iKohl et al.1 119981: 
lAntonucci et al.l 120001) . A similar temperature anisotropy is 
measured in situ at r > 0.3 AU for protons in fast-solar-wind 
streams with p ^ 1, despite the fact that (double) adiabatic 
expansion acts to decrease 7±/7|| dMarsch et ai.ll 19821 l2004i: 

iHelhnger et al.ll2006h . where p = Mp/B'^ is the ratio of the 
plasma pressure to the magnetic pressure. Thus, in coronal 
holes and fast wind with p ^ 1, ions receive > 1/2 of the 
total heating, and ion heating is mostly "perpendicular to the 
magnetic field." 

Because the rms amplitude of the magnetic-field fluctua- 
tion bB at the dissipation scale is <C Bq, the damping of turbu- 
lent fluctuations can be treated, to a first approximation, using 
the Vlasov-Maxwell theory of linear waves. In this theory, 
Alfven waves are virtually undamped when k±Pp <ti 1 and 
CO ^ Hp, where Pp is the rms proton gyroradius and i2p is the 
proton cyclotron frequency. However, as k±Pp increases to 
values > 1, the Alfven waves (AWs) become kinetic Alfven 
waves (KAWs), the ions begin to decouple from the electrons, 
and the waves develop fluctuating electric-field and magnetic- 
field components parallel to Bq (Hasegawa and Chen 197^ 
ISchekochihin et all 120091) . For KAWs with k^ < k^ and 
CO ^ Hp, the primary damping mechanisms are Landau damp- 
ing and/or transit-time damping, which lead to parallel heat- 
ing of the plasma, no t perp endicular heating (Ouataert'1998: 
Leamon et al.l 119991 [Cranmer and van Ballegooijen 200^ 
Garv an d Nishimura 2004). Moreover, in low-p plasmas, the 
waves damp almost entirely on the electrons, because ther- 
mal ions are too slo w to satisfy the Landau resona nce con- 
dition co-/t||V|| = (Oua taert|[l99l lGruzinovl[T998l) . Thus, 
if KAWs damp according to linear Vlasov theory, then they 
are unable to explain the strong perpendicular ion heating that 
is inferred from observations. This discrepancy casts doubt 
on the viability of low-frequency AW/KAW turbulence as a 
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mechanism for heating coronal holes and the fast solar wind. 

A number of studies have gone beyond linear Vlasov 
theory to investigate the possibility of perpendicular 
ion heating by l ow-frequency AW/KAW turbulence. 
iJohnson and Chen^ (lOOf), 'Chen et al." ("lOOB, iWhite et al.1 
(I2002D, and Voitenko and Goossens (2004) investigated the 
dissipation of mono-chromatic KAWs and AWs with CO < fip, 
finding that such waves cause perpendicular io n heating if the 
wave a mpli tude exceed s a minimum threshold. iDmitruk et al.l 
(12004') and lLehe etall (12009) simulated test particles propa- 
gating in the electric and magnetic fields resulting from direct 
numerical simulations of magnetohydrodynamic (MHD) 
turbulence at 0.1 < p < 10. They both f ound perpendicular 
ion heating under some conditions, but iLehe et all ( 12009b 
argued that the perpendicular heating seen in both studies is 
due to cyclotron resonance and does not apply to the solar 
wind because it is an artifact of limited numerical resolution. 
iParashar et al.l ('2009) found perpendicular ion heating in 
two-dimensional hybrid simulations of a turbulent plasma, in 
which ions are treated as p articles and electrons are treated 
as a fluid. In addit i on, iMarkovskii and HoUwed (l2002h 
and IMarko vskii et al. ( 2006') investigated high-frequency 
secondary instabilities that are generated by KAWs near the 
gyroradius scale, and argued that such instabilities may be 
able to explain the observed perpendicular ion heating. 

In this paper, we continue this general line of inquiry and 
address an important open problem: determining the perpen- 
dicular ion heating rate in anisotropic, low-frequency (co < 
ilp), AW/KAW turbulence as a function of the amplitude 
of the turbulent fluctuations at the gyroradius scale. In sec- 
tion |2] we develop a phenomenological theory of stochastic 
ion heating, obtaining an approximate analytic expression for 
the heating rates of different ion species. We also present sim- 
ulations of test particles propagating in a spectrum of AWs 
and KAWs to test our phenomenological theory and to deter- 
mine the two dimensionless constants that appear in our ex- 
pression for the heating rate. In section|3]we apply our results 
to perpendicular ion heating in coronal holes and the fast solar 
wind. 

2. STOCHASTIC ION HEATING BY ALFVENIC TURBULENCE AT 
THE GYRORADIUS SCALE 

We consider ion heating by fluctuations with transverse 
length scales X± (measured perpendicular to Bq) of order the 
ion gyroradius p = v±/Q. (i.e., k±p ~ 1), where Q. = qBo/mc 
is the ion cyclotron frequency, and q and m are the ion charge 
and mass. We assume that p > Pp, where 



(1) 



is the rms proton gyroradius in the background magnetic field. 



l2kBT-p 



(2) 



is the rms perpendicular velocity of protons, Tp is the (per- 
pendicular) proton temperature, and mp is the proton mass. If 
p ^ Pp, then the gyro-scale fluctuations are AWs. If p ~ Pp, 
then the gyro-scale fluctuations are KAWs. 

We define 5vp and dBp to be the rms amplitudes of the fluc- 
tuating velocity and magnetic-field vectors at k±p ^ 1. Sim- 
ilarly, 5Ep and 5<t>p are the rms amplitudes of the fluctuating 
electric field and electrostatic potential at ^^P 1- We as- 
sume that 5vp, 5Bp, 5Ep, and 5<t>p are related to one another in 



the same way that the magnitudes of the fluctuating velocity, 
magnetic field, electric field, and electrostatic potential are re- 
lated in a linear (kinetic) Alfven wave. Thus, since k±Pp<l, 



54>p ^ pdEp, and 



5£„ 



g'84>p 



5vpBo 



5vn 



(3) 



(4) 



The fractional change in an ion's perpendicular kinetic en- 
ergy mv^/2 during a single gyro-period is then given by 



2^5<I>p 



'2e, 



where 



v± ' 



(5) 



(6) 



When e ^ 1, an ion's kinetic energy is nearly constant 
during a single gyro-period. If in addition 5Bp ^ Bq, then 
the ion's orbit in the plane perpendicular to Bq closely ap- 
proximates a closed circle in some suitably chosen reference 
frame. In this case, the ion possesses an adiabatic invariant 
of the form J ~ § pdq that is conserved to a high degree of 
accuracy, where q is the angular coordinate corresponding 
to the particle's nearly periodic cyclo tron gyration a nd p is 
the canonically conjugate momentum (lKruskallll962h . In the 
limit of small e, J is approximately equal to the magnetic mo- 
ment jj — mv\/2B. The near conservation of J implies that 
perpendicular ion heating is extremely weak. In Appendix lAl 
we present a calculation for electrostatic waves with kx_p ^ 1 
and e <C 1 that illustrates how the leading-order terms in the 
time derivative of v\ are unable to cause secular growth in Tx . 

On the other hand, as e increases from to 1, the fractional 
change in an ion's perpendicular kinetic energy during a sin- 
gle gyro-period grows to a value of order unity. We treat the 
spatial variations in the electrostatic potential <I> at fcxP 1 as 
random or disordered, as is the case in turbulence or a spec- 
trum of many randomly phased waves. Thus, when e exceeds 
some threshold (whose value we investigate below), an ion's 
orbit in the plane perpendicular to Bo becomes chaotic. In 
this case, the ion's orbit does not satisfy the criteria for the ap- 
proximate conservation of J (Kruskal 1962), and perpendic- 
ular ion heating becomes possible (Johnson and Chenal2001l : 



iChen et al.ll200U IWhite et al.ll2002 

To estimate the rate at which ions are heated, we begin 
by considering the Hamiltonian of a particle of charge q and 
mass m. 



H = q^+-Lfp^lA 

Zm \ c 



(7) 



where A is the vector potential, and p is the canonical mo- 
mentum. Hamilton's equations imply that 



dH _ 34> q\ dA 
dt ^ 3f c 3? ' 



(8) 



where v — m^' (p — qA./c) is the particle's velocity. The elec- 
tric field is given by E = — V<I> — c^^dA/dt. The second term 
in equation ^ is Es, where Es ~ — c^'3A/3r is the part of 
the electric field that has a nonzero curl. In AWs and KAWs 
with CO < Hp and k±Pp < 1, Eg is negligible compared to 
the total electric field in low-p plasmas [see equation (46) of 
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iHollwed ([T999l) 1. which are our primary focus, and so from 
here on we neglect the second term in equation (O. 

When d(S>/dt > 0, a particle can gain potential energy, ki- 
netic energy, or both. For example, if an ion interacts with an 
electrostatic wave with wavelength ^ p and frequency <C ii, 
then the ion's guiding center drifts with velocity cE x Bq/Bq. 
The particle's kinetic energy undergoes small-amplitude os- 
cillations due to its gyro-motion. However, because its guid- 
ing center moves perpendicular to VO, there is no significant 
secular change in its kinetic energy. The ion's magnetic mo- 
ment /J is almost exactly conserved, and the change in its total- 
energy is almost exactly equal to the change in its potential 
energy. 

On the other hand, if a particle enters a region in which 
d<P/dt > and then leaves this region, moving up and down 
the potential gradient, then it can gain kinetic energy as il- 
lustrated in figure [T] The "wire-mesh" surface in the up- 
per panel of this figure represents <P{xi,X2,t) at some initial 
time, and the lower panel shows <i>{x\,X2,t) at a later time. 
We take the maximum of <I> to be located at a = 0, where 



o = y/xj +X2- We have assumed that 3<I>/3f > at a < Oo 

and d<P/dt = at a > Oo, where Oq is the approximate ra- 
dius in the xi —X2 plane of the "potential-energy hills" that 
appear in the figure. The thick solid line shows the value of 
<t> along the trajectory of a particle moving in a straight line 
in the xi —X2 plane. Because d'P/dt > 0, the potential-energy 
hill is shorter when the particle is "climbing up" and higher 
when the particle is "rolling down." The particle thus experi- 
ences a net gain of kinetic energy from "rolling over the hill." 





Fig. 1. — The potential <I> in tlie x\ —X2 plane at some initial time (upper 
surface), at a later time (lower surface), and along the trajectory of a particle 
moving in a straight line in the jci — X2 plane (thick solid line). The "potential 
energy hill" is shorter when the particle rolls to the top, and higher when the 
particle rolls down, so the particle gains kinetic energy as it rolls over the hill. 

We now estimate the rate at which ions are heated by AW 
or KAW fluctuations with X± ^ p. We note that the condition 
A-x ^ p is intended to encompass structures with X± ~ 0.5p, 
which we invoke below when discussing equation (l23T l. How- 



ever, we ignore fluctuations with X±^ p otX± through- 
out this discussion. Although we are interested in stochastic 
ion orbits, we can still define an effective guiding-center po- 
sition, 

V X b 

R = r+^, (9) 

where b ~B/B and r is the ion's instantaneous position. 
When e ^ 1, the particle gyrates smoothly about position R. 
As e increases towards 1, the particle's motion becomes more 
complicated, but the particle remains within a distance ^ p of 
position R. Taking the time derivative of equation ^ and us- 
ing the equation d\/dt ^ {q/m)(E + \ x B/c), we obtain the 
equation 

dR £ cExB 

:yiib+ +..., (10) 



dt 



B2 



where the ellipsis (. . .) represents terms proportional to deriva- 
tives of B, which we ignore in our approximate treatment. 
During a single cyclotron period, an ion passes through a 
small number of uncorrected fluctuations or "structures" of 
transverse scale ^ p. Within different structures, the vector 
cE X B/B^ has a similar magnitude 5vp) but points in dif- 
ferent directions. The time average of cE x B/B^ over a sin- 
gle cyclotron period is thus somewhat smaller than, but of 
order, Svp. The time At required for an ion's guiding center to 
move a distance p is thus approximately 



At' 



_P_ 

Sv'o 



(11) 



[In writing equation (fTTI ). we have assumed that the gyro-scale 
fluctuations do not oscillate on a time scale ^ At, and we con- 
tinue to make this assumption in the analysis to follow.] Each 
time the particle moves a distance p perpendicular to Bo, it 
encounters different and uncorrected gyro-scale electromag- 
netic fields. Thus, dR/dt decorrelates after a time Af, and the 
particle's guiding center undergoes a random walk in space 
with diffusion coefficient ^ p^/At. 

Similarly, when £ is sufficiently large that the ion's motion 
becomes stochastic, the value of dH/dt decorrelates after a 
time At, and the particle undergoes a random walk in energy. 
In contrast, as shown in AppendixlAl as £ — > the interaction 
between ions and gyro-scale electrostatic -potential structures 
is not a Markov process; instead, changes in H are correlated 
over long times, and to leading order in £ are reversible and 
bounded. Returning to the stochastic case, we define 3<I>/3f 
to be the rms value of 3<I>/3f associated with fluctuations with 
~ p. The rms change in H during a time At is then 



a* 

AH r^q — At. 
dt 



(12) 



An ion undergoing stochastic motion can gain kinetic energy 
in the same way as the particle illustrated in figure [T] If the 
ion spends a time Af localized within a flux tube of cross- 
sectional area ^ p^ and length ^ |v|| \At, it exits this flux tube 

in a random direction. Thus, if 3<i>/3f is on average posi- 
tive during this time interval within the flux tube, it does not 
follow that the ion will move to a region of larger <!> after a 
time At, where <i> is the electrostatic potential associated with 
fluctuations with A-^ ~ p. On the contrary, the change in <!> 
along the ion's path is only loosely correlated with the aver- 
age change in <!> within the flux tube. As a result, the change 
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in the ion's kinetic energy during a time A? is of the same 
order of magnitude as the change in its total energy given in 
equation (fT2b.^ Because V<P is nearly perpendicular to B, and 
because the ion's guiding center moves perpendicular to B by 
a distance of order A-^ ~ p during a time At, the ion's perpen- 
dicular kinetic energy K±^ = mv^ / 2 changes by an amount of 
order 

AK^r^AH (13) 

during a time Af. We discuss the parallel kinetic energy 
following equation ( l24l l below. We define an effective fre- 
quency COeff for gyro-scale fluctuations through the equation 



(14) 



For example, if the gyro-scale fluctuations consist of waves 
with a single frequency CO, then cOeff = CO. With the use of 
equations (|4|l and (fTTT l. we can rewrite equation ( fT3] l as 



AKi 



L COeff p. 



(15) 



The kinetic-energy diffusion coefficient Dk ^ {AK±^^/At is 
then given by 

Da: -m^vico^ffSvpp. (16) 

When a single ion undergoes kinetic-energy diffusion, the 
ion has an equal likelihood of gaining or losing kinetic en- 
ergy during each "random-walk step" of duration At. On the 
other hand, if a large population of ions undergoes kinetic- 
energy diffusion, and if the phase-space density / of ions is 
a monotonically decreasing function of K±, then the average 
value of K± increases steadily in time. To distinguish between 
properties of individual particles and rms quantities within a 
distribution, we define v^i to be the rms perpendicular ve- 
locity of the ions, which is related to the perpendicular ion 
temperature T± by the equation 



2^b71 



(17) 



(18) 



We also define the rms ion gyroradius. 

We define 5vi to be the rms amplitude of the fluctuating fluid 
velocity at ^ Pi, and we set 

Ei = — . (19) 

For protons, we define 5vp (5Bp) to be the fluctuating fluid 
velocity (magnetic field) at X± ~ Pp, and we define 



6Vn 



(20) 



The time scale for the average value of K± in a distribution of 
ions to double is then roughly 



t\ 



±i 



(21) 



' In contrast, in the small-e limit addressed in Appendix |A] the change 
in a particle's total energy is almost exactly equal to the change in the gyro- 
averaged potential energy. 



where Dk{ is the value of Dk for ions with = vxi and p = 
Pi. The perpendicular ion heating rate per unit mass is then 



■Jh, or 



co; 



'eff,i 



SviPi, 



(22) 



where cOgff.i is the value of cOeff at p = pi. 

We now consider what determines the value of cOeff in 
anisotropic AW or KAW turbulence. If the turbulence is 
driven at an "outer scale" Lq that is 3> p, the advection or 
"sweeping" of structures with A-^ ^ p by the outer-scale ve- 
locity fluctuations leads to rapid time variations in <I> at a fixed 
point in space. On the other hand, these large-scale velocity 
fluctuations advect both the ions and the small-scale struc- 
tures in the electric and magnetic fields. Thus, if one consid- 
ers ions within a flux tube of radius ^ p and length ^ Lq, 
and if one transforms to a frame of reference moving with the 
average velocity of that flux tube, then the rapid time varia- 
tions resulting from large-scale advection disappear This in- 
dicates that large-scale sweeping does not control the rate of 
ion heating or the value of cOeff in equation (|22] |. On the other 
hand, electrostatic-potential structures at scale A,^ ~ 0.5p are 
advected by velocity fluctuations at the same scale, and there 
is no frame of reference in which the velocities at Ax — 0.5p 
vanish at all points along an ion's gyro-orbit. This advection 
by velocity fluctuations with Ax — 0.5p causes 3<I>/3f to have 
a value of ^ 5<t>p5vp/p, which gives 



COeff ' 



8vd 



(23) 



where we have neglected factors of order unity, such as the 
ratio between Svp and the rms amplitude of the velocity 
fluctuation at Ax — 0.5p. Put another way, the advection 
of electrostatic-potential structures at Ax ~ 0.5p, which are 
rooted in the electron fluid, leads to a partial time derivative 
of <I> that ions can feel, and which energizes ions through the 
process illustrated in figure [T] We note that in "imbalanced" 
(or cross-helical) AW turbulence, in which the majority of the 
waves propagate either parallel to Bo or anti-parallel to Bo, 
the energy cascade time for the majority waves can greatly 
exceed CO^j^', since the majority waves are cascaded by the 
smaller-amplitude waves propagating in the opposite direc- 
tion. Nevertheless, even for imbalanced turbulence, the argu- 
ments leading to equation ( l23T l continue to hold. 

As discussed following equation (|6]l and in Appendix lAl 
when e is sufficiently small, the changes in H remain cor- 
related (and largely reversible) over long times, so that the 
perpendicular heating rate is strongly reduced relative to our 
estimate in equation (|22] |. To account for this, we introduce 
a multiplicative suppression factor onto the right-hand side 
of equation (|22] | of the form exp(— C2/ei). We also add an 
overall coefficient ci to the right-hand side of equation (l22l l to 
account for the various approximations we have made. Both 
c\ and C2 are dimensionless constants whose values depend 
upon the nature of the fluctuations (e.g., whether the fluctua- 
tions are waves or turbulence, the type of turbulence, etc) and 
the shape of the ion velocity distribution. Substituting equa- 
tion ( l23T l into equation (l22T i. we obtain 



y± — — exp 



(24) 



We emphasize that for protons in low-p plasmas, Ep ~ 
P^'/^SBp/Bo, and thus ep can approach unity even if dBp/Bo 
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remains <^ 1 . 

The change in an ion's parallel kinetic energy K^^ = mv|/2 
during a time Af due to the parallel electric field E^^ is AK^^ ^ 
qE\\vi^\At. We have restricted our analysis to AWs and KAWs 

with A,^ ^ p ?J Pp and co^' > Af ^ A-^/5vp. This condition on 
the wave frequency implies that the parallel wavelengths of 
such fluctuations satisfy the inequality > pv^/Svp A-^- 

When Me/mp < p < 1 and A,^ > Pp, E\\/E±_ ^ Pp/(AxA||), 
where is the electron mass dHoUwegl 119991) and E± is 
the electric-field component perpendicular to B. Thus AK\\ is 
< V|| /va times the value of AK± in equation (fTST l. For thermal 
ions in low-P plasmas, V|| ^ va- Thus, when Ei is sufficiently 
large that stochastic heating is important, stochastic heating 
leads primarily to perpendicular ion heating rather than par- 
allel heating. For AWs/KAWs in low-p plasmas, the parallel 
component of the magnetic mirror force is much less than qE\\ 
( lHollweglll999h and thus does not affect our conclusions re- 
garding anisotropic heating at p <C 1 . 

2.1. Test-Particle Simulations of Proton Heating 

To test the above ideas, we have numerically simulated 
test-particle protons interacting with a spectrum of randomly 
phased KAWs. The protons' initial locations are chosen ran- 
domly from a uniform distribution within a volume encom- 
passing many wavelengths perpendicular and parallel to Bq. 
The protons' initial velocities are drawn randomly from an 
isotropic Maxwellian distribution of temperature Tp. For each 
particle, we solve the equations 



equation ( l23T l for KAW turbulence at ^^Pp ^ 1- We then set 



dx 
dt 



and 



d\ 
'di 



q 



E- 



vxB 



(25) 



(26) 



using the Bulirsch-Stoer method (' Press et anil992 >). We take 
B = Bqz + Bi, where Bq is constant. We take E and Bi to be 
the sum of the electric and magnetic fields from 162 waves 
with randomly chosen initial phases, with two waves at each 
of 81 different wave vectors. At each wave vector, there is 
one wave with Co/A;^ > and a second wave with (o/k- < 0. 
This second wave has the same amplitude as the first, so that 
there are equal fluxes of waves propagating in the +z and ~z 
directions. The 81 different wave vectors consist of 9 wave 
vectors at each of nine different values of k^, denoted k^j. 
The k^j can be expressed in terms Pp. In particular, the val- 
ues — ln(A:x;Pp) are uniformly spaced between —4/3 and 
4/3; i.e., \|/; = -4/3 + j/3, with ;= 0, 1, ... ,8. We regard 
the values as corresponding to cell centers in a uniform 
grid in \|/ = In(fcxPp)^ with grid spacing Ax|/ = 1 /3. The mid- 
dle three grid cells, with j — 3, 4, and 5, thus correspond to an 
interval of width unity in ln(A:x) space centered on A:xPp = 1. 
We define the rms amplitudes of the gyro-scale velocity and 
magnetic-field fluctuations 5vp and 5Bp in our simulations by 
taking the rms values of the E x B velocity and magnetic-field 
fluctuation resulting from the KAWs in these middle three 
grid cells. At each k^j we include 9 different values of the az- 
imuthal angle i^'mk space, (|)/ — 27iZ/9, where / = 0, 1 , . . . , 8. 
At each k^j there is only a single value of which we de- 



note k^i^j. We choose k^ so that the frequency at k±^ = k^/[ and 
^11 ~ ^||4 equals ki^^hvp. The linear frequency of our gyro- 
scale KAWs is thus comparable to the value of COeff given in 



21/ 
ku 



{k^j/k^^fl^ if0<j<4 
{k^j/k^^yl^ if4<;<8 



(27) 



Our formula for k\\j at y < 4 is chosen so that the wave pe- 
riods are comparable to the energy cascade time s c ales i n 
the critical-balance theory of iGoldreich and Sridhaii (119951) . 
while the formula for j > 4 is chosen so that the wave periods 
match the nonlinear time scales in the critical-balance theory 
of ICho and LazarianI (2004). All waves at the same fc^ have 
the same amplitude, and (since there are the same number of 
waves at each k^j) we take the amplitude of the magnetic- 

-1/3 

field fluctuation in each wave to ht °^ k^' for A:xpp < 1 

-2/3 

and 0= A: I ' for A:xPp > 1, again motivated by the theories of 
Goldreich and Sridhar (1995) and Cho and Lazarian (2004). 

The relative amplitudes of the different components of E 
an d Bi for each wave are taken from the two-fluid theory 
of iHollwe^ (Il999h . To apply this theory, we choose plasma 
parameters that are characteristic of coronal holes. In par- 
ticular, we set Pe = 8%nkBTe/Bl = 0.003, va = 0.003c, and 
Te — 0.5rp, where n is the electron number density (equal to 
the proton number density), and Te is the electron tempera- 
ture. 

Using the above procedures, we have carried out seven sim- 
ulations with different values for the overall normalization of 
the wave amplitudes, with 5Bp/Bo ranging from 4.8 x 10^^ 
to 1.9 X 10^^. Given the polarization properties of KAWs 
and our method for constructing the wave spectra, the value 
of 5vp/vA is 1.19 times the value of dBp/Bo in each simula- 
tion. The wave frequencies reach their maximum values in the 
largest-5Bp/Bo simulation. In this simulation, co = 0.29£2p at 
A:xPp = 1, and co = 0.82i2p at the maximum value of k±Pp, 
which is 3.79. Although this maximum frequency is close 
to i2p, the cyclotron resonance condition CO — A;||V|| — IClp 
(where / is any integer) is not satisfied, because the parallel 
thermal speed of the protons is only 0.055va and ^ k±. 
For most of the waves in these simulations, co <C iip. 

We determine the perpendicular proton heating rate per unit 
mass 2±p in the simulations by plotting (v^) versus time, fit- 
ting this plot to a straight line to determine {d/dt){v^), and 
then setting 2±p = 0.5{d/dt){v^j), where (...) indicates an 
average over the 10^ particles in each simulation. When fitting 
the plot of (v^) versus time, we ignore the first 10 cyclotron 
periods, because during the first couple gyro-periods the par- 
ticles undergo a modest apparent heating as they "pick up" 
some portion of the E x B velocity of the waves. We find that 
after (v^) increases by between 20% and 40%, the heating 
rate starts to decrease for two reasons. First, the small-v'x part 
of the velocity distribution flattens, after which this part of the 
distribution is no longer heated as effectively. Second, as (v^) 
increases, £p decreases. We neglect this later stage of weaker 
heating when constructing our fits to the (v^lO) plots, so that 
the measured heating rates correspond to Maxwellian distri- 
butions. (For the smallest values of 5vp, we do not observe a 
second stage of weaker heating, because the test-particle ve- 
locity distributions do not change very much during the sim- 
ulations, which last IQ'^Q.p^.) We illustrate this procedure in 
figure |2] for a run with 5vp/vxp = 0.15. In this case, we de- 
termine Q±p from the slope of the long-dashed line, which is 
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2000 



Fig. 2. — i'^ and i'^ averaged over the 10 particles in a simulation with 

Si'p/v^p = 0.15, = 0.003, va = 0.003c, and = O.STj,. The two solid- 
line curves correspond to our basic numerical method. We determine Qxp in 
this simulation from the slope of the long-dashed line. The short-dashed fine 
shows (vy ) in a modified simulation with the same parameters in which E is 

replaced by E' = E + b(£:j - b ■ E) . 



our fit to the (v^) data over the interval 10 < Q.pt < lO"'. 

In figure |3] we plot the values of Q±p for several different 
values of Ep. Each x in this figure corresponds to a separate 
simulation with a different value of 8vp but the same initial 
proton temperature. The solid line is the proton heating rate 
from equation ( l24l i with ci = 0.75 and C2 = 0.34; that is, 



0.75 (8v 



Pp 



exp 



0.34 



(28) 



We expect the constants ci and C2 to be fairly insensitive to 
variations in pe, 7p/7e, and va/c (at least within the range 
of solar-wind-relevant parameters), in which case Q^^ de- 
pends on the plasma parameters primarily through the ex- 
plicit Pp and Ep terms in equation (|28] |. The values of ci 
and C2 in equation ( |28] | presuppose the presence of a broad 
spectrum of AWs and KAWs bracketing the perpendicular 
wavenumber k± — (Pp)^', encompassing at a minimum the 
range 0.3 < k±Pp < 3. A spectrum of at least this width is 
probably present in the solar wind, the only uncertainty being 
the value of the dissipation wavenumber beyond which the 
wave power spectrum decreases exponentially with increas- 
ing k±. If the simulations described in this section are re- 
peated without the smallest three values of k± and without 
the largest three values of k± (keeping the wave amplitudes 
fixed at the middle three values of k±), then the proton orbits 
become less stochastic, and g±p decreases significantly. (The 
exact amount by which 2±p decreases depends upon the value 
of 8p.) We have omitted waves at k±Pp < 0.26 and k±Pp > 3.8, 
but we expect that waves at such scales do not have a strong 
effect on perpendicular ion heating, provided co is sufficiently 
small that the cyclotron resonance condition can not be satis- 
fied. It is possible that in some cases strongly turbulent fluctu- 
ations with k±Pp ^ 1 and nonhnear time scales ~ i2p ' could 
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Fig. 3. — Numerical results (plotted with xs) for the perpendicular heat- 
ing rate Q^p for test-particle protons interacting with a spectrum of randomly 
phas ed K AWs. The solid line is equation (28), and the dashed line is equa- 
tion j28t with the "^u conservation" factor exp(— 0.34/ep) replaced with unity. 



heat ions through a broadened cyclotron resonance, but a de- 
tailed investigation of this process is beyond the scope of this 
study. 

We reiterate that the values of ci and C2 in equation ( l28T l are 
not universal, but instead depend on the type of fluctuations 
that are present. In true turbulence (as opposed to random- 
phased waves), the value of C2 may be smaller than in our 
simulations (indicating stronger heating), because a signifi- 
cant fraction of the cascade power may be dissipated in co- 
herent structures in which the fl uctuating fields are larger than 
their rms values (D mitruk et alj 2004). 

The lower solid-line curve in figure|2]plots (v|) versus time 

in the simulation with ep = 0.15, = 0.003, va = 0.003c, 
and Te = 0.5Tp. During the interval 10 < £2pf < 2200, the in- 
crease in (v|) is about one-fourth the increase in (v^^). How- 
ever, most of the increase in (v|) is an artifact of our numerical 

method, which equates the parallel electric fields of the waves 
with the z component of the electric field in the simulation, 
and the perpendicular electric field of the waves with the x and 
y components of the electric field in the simulation. The local 
magnetic field in our simulations, however, is not parallel to 
the z axis, but instead has nonzero x and y components result- 
ing from the magnetic-field fluctuations. As a result, part of 
the perpendicular wave electric field is converted into a par- 
allel electric field in the simulation, artificially enhancing the 
parallel electric field seen by the particles. To eliminate this 
effect, we have repeated this simulation replacing the local 
electric field E seen by each particle with the adjusted elec- 
tric field E' = E + b(£', - b • E), where b B/B and B is 
the local value of the magnetic field. In this new simulation, 
the parallel electric field b • E' is the sum of the parallel elec- 
tric fields of the individual waves in the simulation and does 
not include any contribution from the perpendicular electric 
fields of the individual waves. The value of (v|) in this mod- 
ified simulation, shown as a dashed fine in figure |2l does not 
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increase significantly during the course of the simulation (in 
fact it decreases slightly), consistent with our argument above 
that parallel heating is weak when p <C 1 . 

2.2. Proton Heating at ki^p^ \ as a Fraction of the 
Turbulent Cascade Power 

The cascade power per unit mass at fc^pp ^ 1, which we 
denote Y, depends upon whether the turbulence is "balanced" 
or "imbalanced," where balanced (imbalanced) turbulence in- 
volves equal (unequal) fluxes of waves propagating parallel 
to Bo and anti-paraUel to Bq. In balanced KAW turbulence, 



8Bp 
Bo 



(29) 



where Ck is a dimensionless constant (Howes et al. 2008a). It 
can be inferred from the numerical simulations of Howes et al. 
(2008b) that Ck = 2.0 (G. Howes, private communication). In 
the simulations of section [ZTl 5Bp/Bo = 0.845vp/vA, and we 
make the approximation that this same ratio is characteristic 
of KAW turbulence in general. Combining equations ( |28] ) and 
( |29l ). we obtain 



e±p 

r 



3.0 exp 



0.34 



(30) 



We expect that Ck, like the constants c\ and C2, depends only 
weakly on P, Tp/Te, and va/c (at least for solar- wind-relevant 
parameters), so that the numerical constants 3.0 and 3.4 in 
equation ( l30l l are relatively insensitive to the plasma parame- 
ters. Equation ( l30l l implies that perpendicular proton heating 
by KAWs with ^^Pp 1 absorbs > 1/2 of the cascade power 
at k±Pp 1 when Ep exceeds 

ecrit = 0.19. (31) 

The cascade power in imbalanced AW turbulence is smaller 
than in balanced AW turbulence with the same total fluctu- 
ation energy, because the AW energy cascade requ ires in- 
teractions between opp ositely propagating waves (Iroshnikov) 
IT96llKraichnanllT96l . At k^Pp - 1, KAWs propagating in 
the same direction can interact nonlinearly with one another, 
but the importance of such interactions relative to interactions 
between oppositely propagating waves is not well known. De- 
spite this uncertainty, we expect that if AW/KAW turbulence 
is imbalanced at fe^Pp 1> then the cascade power at ^^Pp 1 
is less than in equation ( |29] l. On the other hand, it is unlikely 
that imbalance strongly affects Q±p if 5vp is held fixed (ex- 
cept for particles with V|| ~ ±va, as discussed in section |24] |. 
We thus expect perpendicular proton heating to absorb at least 
50% of the cascade power at k±Pp ~ 1 in imbalanced turbu- 
lence even when ep is somewhat smaller than 0.19. 

2.3. Proton Heating versus Electron Heating by KAWs with 

k^Pp - 1 

Stochastic proton heating removes energy from KAW fluc- 
tuations with kx_Pp ^ 1, resulting in an effective damping rate 
for these fluctuations, which we denote Yp. The value of Yp is 
given by the relation 

2YpEw = exp, (32) 

where Ew is the energy per unit mass of the KAW fluctua- 
tions at ^^Pp ^ 1. The factor of 2 in equation (l32T i is included 
to make Yp analogous to a linear wave damping rate, in the 
sense that the rate at which linear waves lose energy is twice 



the product of the damping rate and the wave energy. To esti- 
mate the value of Yp in AW/K AW turbulence, we use the test- 
particle calculations in section|2T|for a spectrum of randomly 
phased KAWs. We take to be the energy per unit mass of 
the full spectrum of waves in these simulations. (This choice 
leads to a conservative estimate of Yp, since the damping is 
likely concentrated in the subset of the waves with k^fip > 1 .) 
On the other hand, we continue to define (5vp)^ as the mean- 
square E X B velocity associated with KAWs with values of 
kx^ lying within a logarithmic interval of width unity centered 
on ^^Pp = 1. With these definitions, Ew = 2.1 (5vp)^ in all of 
the simulations in section IZTI Combining equations ( |28] ) and 
( |32] |. we obtain 



Yp = 0.18epnp exp 



0.34 



(33) 



In low-P plasmas, small-amplitude KAWs with k^^Pp — 1 
and CO <C iip undergo electr on Landau damping but negligi- 
ble linear proton damping (lOuataerti 1 1 998tlGruzinovl IT 9981 : 
iGarv and Nishimurall2004l) . Using the numerical method de- 
scribed by'Ouataert (1998j and Howes et al. (2008a), we nu- 
merically solve the full hot-plasma dispersion relation to find 
the electron damping rate Ye of KAWs with k^Pp = 1 and 
CO <C £ip for a range of values of k\\, Tp/T^, pp, and va/c, where 
Pp = SnnkBTp/Bl. We find that if me/ nip < pe < 1, va < c, 
and 0.1 < Tp/Te < 10, then the damping rate at k±Pp = 1 is 



well fitby the formula Ye =9.5 x 10"^(7;/rp)i/2pp 
or equivalently 



-1/2 



Ye = 9.5 X 10-\x' 



PpTp 



1/2 



^IIVaI, 



(34) 



where % = fc^Svp/|A:| |VA|- In some theorie s of strong MHD 
turbulence % ^ I ( Goldreich and Sridha^ 119951: iBoldvrevI 
I2OO6) . This condition, some times referred to as critical 
balance, may characterize AW/KAW fluctuations in coronal 
holes and the solar wind at k±pp ^ 1. On the other hand, if 
the frequencies of the waves launched by photospheric mo- 
tions are sufficiently small, then AW/KAW turbulence at a 
heliocentric distance of a few solar radii may be more "two- 
dimensional" than in critical-balance models, with smaller 
values of ^|| and a larger value of %. 
Combining equations ( |33] ) and (l34l i. we obtain 



PpTp 



1/2 



exp 



0.34 



(35) 



The ratio Yp/Ye approximates the ratio of the proton heating 
rate to the electron heating rate resulting from KAW fluctua- 
tions at k±pp ~ 1 in the low-p conditions present in coronal 
holes and the near-Sun solar wind. (At P > 1, linear KAW 
damping on the protons becomes important, increasing the 
proton heating rate.) We note that if the damping time scales 
Yp ^ and Ye^' are both much longer than the energy cascade 
time at k±Pp ^ 1, then most of the fluctuation energy will cas- 
cade past the proton-gyroradius scale to smaller scales. In that 
case, the division of the turbulent heating between protons and 
electrons will depend primarily upon how fluctuations dissi- 
pate at ^xPp 3> 1. 

2.4. How the Heating Rate Depends on q, m, P, and vh/va 
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If we re-run our simulations, keeping only waves with 
Co/A:^. > 0, and consider a thermal distribution of test-particle 
protons with a nonzero average velocity equal to vaZ, then 
the perpendicular heating rate is strongly reduced. This is 
because the electric field of an Alfven wave (or KAW with 
X± ~ Pp) vanishes (or is strongly reduced) in a reference 
frame moving at speed va in the same direction as the wave 
along the background magnetic field. This effect may explain 
the observation that the perpendicular heating of a particles 
in the solar wind is reduced when the differential flow ve- 
locity of a particles relative to protons (in the anti-Sunward 
direction) approaches va (Ka sper et al.l i2008). at least in re- 
gions where anti-Sunward propagating KAWs dominate over 
Sunward-propagating KAWs. 

If we hold 5Bp/Bo fixed but increase pp to 1, then the per- 
pendicular proton heating rate is dramatically reduced, be- 
cause Ep = 5vp/v± ~ P^'/^5Bp/Bo decreases by a large factor. 
On the other hand, the protons in these Pp 1 simulations un- 
dergo signific ant parallel heat ing, consistent with results from 
linear theory ('Ouataert"1998') and recent test-particle simula- 
tions of ions propagating in numerically simulated MHD tur- 
bulence (iLehe et al.ll2009l) . 

If we re-run our simulations but use ions instead of 
protons (but with the same temperature as the protons), then 
the perpendicular heating rate is much larger. This is in large 
part because e is larger for (and other heavy ions) than for 
protons at the same temperature, a point to which we return 
in section |3] Another reason for enhanced heavy-ion heat- 
ing can be seen from equation ( 1211 1. We rewrite this equation 
with the aid of equation ( |23] |. increasing fi by exp(c2/ei) for 
the same reasons that we reduced Q± by this same factor in 
equation (l24l l. to obtain 



viiPi 



exp 



(36) 



In a number of theories of MHD turbulence, the ratio 
(5vpi)-'/pi is relatively (or completely) insensitive to the value 
of Pi, provided Pi is in the inertial range of the turbulence. On 
the other hand, for ion species at equal temperatures, v^j is 
inversely proportional to the ion mass. Thus, even aside from 
the exponential factor in equation (l36l l. the heating time scale 
is shorter for heavier ions than for lighter ions at the same 

temperature. 

Finally, if we repeat the simulations of section IZTI for test- 
particle ions with Pi ^ Pp, and with values of k±Pi centered 
on 1 so that the gyro-scale fluctuations are now AWs, we re- 
cover similar values for the perpendicular heating rate per unit 
mass. Stochastic perpendicular ion heating thus does not re- 
quire the particular polarization properties of KAWs, but op- 
erates for both KAWs and AWs, as we have argued in our 
heuristic derivation of equation (l24l i. 

2.5. Lack of Perpendicular Heating by AWs with k±p[ <^ 1 

In turbulent flows, the rms variation in the velocity across a 
perpendicular scale A.^, denoted dvx^, typically increases as 
some positive power of X± when X± is in the inertial range. 
As a result, the variation in the electrostatic potential across an 
ion's gyro-orbit is dominated by the fluctuations at the large- 
scale end of the inertial range, suggesting that these large- 
scale fluctuations might make an important contribution to 
the perpendicular heating rate. This suggestion, however, is 
incorrect, because AWs with k±p ^ I cause an ion's guid- 
ing center to drift smoothly at velocity cE x B/B^, but do not 



cause an ion's motion to become chaotic. If one transforms to 
a reference frame that moves at the velocity cE x B/B^ evalu- 
ated at the ion's guiding-center position, then the variation in 
q<P across the ion's gyroradius is a small fraction of mv^/2. 
The ion's trajectory in the plane perpendicular to Bq in this 
frame is approximately a closed circ le, and the ion 's magnetic 
moment jj is then nearly conserved (iKruskalll 19621) . 

3. PERPENDICULAR ION HEATING IN CORONAL HOLES AND THE 
FAST SOLAR WIND 

As shown in the previous section, the stochastic ion heat- 
ing rate is a strongly increasing function of Ei = Svi/v^i. For 
fixed turbulence properties, the value of Ei depends upon the 
ion charge q = Ze, the ion mass m = Am^, and the perpen- 
dicular ion temperature T±. For example, if we take the rms 
amplitude of the turbulent velocity fluctuation at transverse 
scale X± to be given by 



= av'A 1 — 



(37) 



for Pp < A-x < Lq, where a and a are dimensionless constants 
and Lq is the outer scale or driving scale of the turbulence, 
then 



Ei = a ( : 



(l-„)/2^(l+«)/2 



(38) 



where dp = v'A/iip is the proton inertial length, Pp — 
SnnpkBTp/B^, n^ is the proton density, Tp is the proton tem- 
perature, and /jp is the mean molecular weight per proton; 
that is, the mass density is /jpnpnip, and the Alfven speed is 

B I ^47i/jp«pmp. If the velocity power spectrum ' is k']^^ 



for Lq ' < < Pi \ then 



-1 



C3 



(39) 



To investigate the possible role of stochastic ion heating 
in coronal holes and the fast solar wind, we evaluate equa- 
tion ( |38] | as a function of heliocentric distance r using ob- 
servationally constrained profiles for the density, temperature, 
and field strength. W e take «p to be given by equation (4) of 
iFeldman et al.l (Il997h . which describes coronal holes out to 
several solar radii, plus an additional component proportional 



to r 



np(r) = 

where d 
We set 



3.23x10** 2.51x10'^ 1.85x10' 



^2 



;15.6 ^3.76 

r//?0. Equation ( l40b gives Wp 

l-(2/3)exp(-t//1.5) 



: 4 cm" 



Tp = 3 x 10 



K- 



(40) 
^ at 1 AU. 



(41) 



(l+0.1t/)<'-8 

which leads to a proton temperature that is 10^ K at the 
coronal base, between 2 x 10^ K and 3 x 10^ K in coronal 
holes, and ~ 2 5 x 10^ K at 1 AU. We take the magnetic field 
strength to be jHollweg and Isenberd2002l) 



l-5(/max-l) 
^6 



1.5 



Gauss, 



(42) 



with /max (the super-radial expansion factor) equal to 5. We 
determine the rms amplitude of the fluctuating wave veloc- 
ity at the outer scale, Svl^ = avA, using the analytical model 
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of IChandran and Hollwed ( 120091) . which describes the prop- 
agation of low-frequency Alfven waves launched outward 
from the Sun, taking into account non-WKB wave reflec- 
tions arising from Alfven-speed gradients as well as the cas- 
cade and dissipation of wave energy arising from nonlinear 
wave-wave interactions. In particular, we set 5vig equal to 
the value of Svi-ms plo tt ed wi th a solid line in figure 6 of 
IChandran and Hollwegl (120091) (the curve corresponding to 
their "extended model"). We take Lq to be 10^ km at the 
coronal base [the limit t/ ~> 1 in equation (l42l) 1. and to be 
proportional to Z?^'/^. 

We consider three different values for the spectral in- 
dex cy. 5/3, 3/2, and 6/5. The value C3 = 5/3 is sug- 
gested by in situ measure ments of magnetic-field fluctu- 
ations in the solar wind (Matth aeus and GoldsteinI 119821: 
iBruno and Carbori3 r2005). as well as some theoretical and 
numerical studies of M HD turbulence (iGoldreich and Sridhaii 
Il995t ICho and Vishniac.,200Q) . The value C3 = 3/2 is mo- 
tiva ted by a different set of theoretical and numerical stud- 
ies dBoldvrevI 120061; iMason et al.ll2006t iPerez and BoldvrevI 
l2009r as well as recent in situ observations of the veloc- 
ity power spectrum in the solar wind (iPodesta et al.l 120071 
iPodesta and Bhattacharjed l2009h . The third value, C3 = 
1 .2, follows from recent numerical simulations of reflection- 
driven Alfven-wave turbulence in coronal holes and the fast 
solar wind (Verdini et al. 2009^?). In these simulations, C3 = 
1 .2 at r < 1 .2Rq, and C3 gradually increases towards 5/3 with 
increasing r. 

In figure m we plot £; for H+, He++, and 0+^ assuming 
fjp = 1.2. Although alpha particles and minor ions are ob- 
served to be hotter than protons in the fast solar wind, we 
have set all the ion temperatures equal to Tp to investigate the 
relative heating rates of different ion species that start out at 
the same temperature. Figure |4] illustrates the general point 
that Ei depends strongly on the spectral index C3. In partic- 
ular, decreasing C3 by 28% from 5/3 to 1.2 increases Ei by 
a factor of > 10 at all radii shown for all three ion species. 
Because Q± depends strongly on £i, Q± is extremely sensi- 
tive to the value of C3. For example, for protons, if C3 = 1.2, 
then Ep ~ 0.2 except at r < l.SR^. The approximations lead- 
ing to equation (|30] | imply that Q±p/r = 0.55 when Ep = 0.2 
(where F is the cascade power at fc^Pp ^1), indicating that 
perpendicular proton heating absorbs a substantial fraction of 
the turbulent heating power when C3 < 1 .2. On the other hand, 
if C3 = 5/3, then Ep < 0.02 and Q±p/r in equation (|30] | is 
< 1.2 X IQ-l 

A second general point illustrated by figured is that when 
C3 is fixed, Ei depends only weakly on r for 2Rq < r < 1 AU. 
As a result, given our assumptions, a large radial variation 
in E within this range of r requires a radial variation in the 
spectral index C3 . As mentione d above, the numerical simula- 
tions of I Verdini etaTI (l2009fcl) found C3 ~ 1.2 at r < I.IR^, 
with C3 increasing towards 5/3 with increasing r. In ad- 
dition, radio observations show that the density-fluctuation 
power spe ctrum is significantly fl atter at r — SRp, than at 
r > l ORp) dMarkovskii and HoUweg 2002; H armon and ColesI 
12005b . These observations and numerical simulations raise 
the possibility that C3 is significantly smaller (and that Q±p/r 
is much larger) close to the Sun than at 1 AU. However, 
the inertial range of reflection-driven AW turbulence in coro- 
nal holes is still not well understood. Likewise, the relation 
between the density power spectrum and the velocity power 
spectrum in the imbalanced AW turbulence found in coronal 



holes is not clear. The r-dependence of C3 thus remains un- 
certain. 

A third point illustrated by figure |4]is that Ei is significantly 
larger for He++ and than for protons at the same temper- 
ature. For example, at equal temperatures, protons and alpha 
particles have the same gyroradius, and Ea = 2Ep, where Ea 
is the value of Ei for alpha particles. Because of the strong 
dependence of Q± on Ei, it is possible that the perpendicular 
heating rate per unit volume from stochastic heating by gyro- 
scale fluctuations is larger for alpha particles than for protons, 
even though Helium comprises only ^ 20% of the mass in 
the solar wind. Depending on the values of 02 and C3, it is 
also possible that Helium absorbs a significant fraction of the 
turbulent cascade power in the solar wind. In addition, the 
comparatively large value of Ei for may explain why 
ions are ob served to be so rnuch hotter than protons in the so- 
lar corona (iKohl et al.lll998l;lAntonucci et al.ll2000h . and like- 
wise for other minor ions. 

4. CONCLUSION 

When an ion interacts with turbulent AWs and/or KAWs, 
and when the amplitudes of the fluctuating electromagnetic 
fields atX± p are sufficiently large, the ion's orbit becomes 
chaotic, and the ion undergoes stochastic perpendicular heat- 
ing. The parameter that has the largest effect on the heating 
rate is E = 5vp/v±, where 5vp is the rms amplitude of the ve- 
locity fluctuation at X± ~ p. In the limit E — >^ 0, the ion's 
magnetic moment is nearly conserved, and perpendicular ion 
heating is extremely weak. On the other hand, as E increases 
towards unity, magnetic moment conservation is violated, and 
stochastic perpendicular heating becomes increasingly strong. 

Using phenomenological arguments, we have derived an 
analytic formula for the perpendicular heating rate Q± for 
different ion species. This formula (equation (l24l i) contains 
two dimensionless constants, ci and C2, whose values de- 
pend on the nature of the fluctuations (e.g., waves versus 
turbulence, the slope of the power spectrum) and the shape 
of the ion velocity distribution. Using test-particle simula- 
tions, we numerically evaluate these constants for the case 
in which a Maxwellian distribution of protons interacts with 
a spectrum of random-phase AWs and KAWs at perpendic- 
ular wavenumbers in the range 0.264 < ^^Pp < 3.79, where 
Pp is the rms proton gyroradius in the background magnetic 
field Bq. The particular form of the wave power spectrum 
that we choose for these si mulations is motivated by the 
"critical b alance" theories o f IGoldreich and Sridharl (119951) 
and ICho and LazarianI ((2004). For this case, ci = 0.75 and 
C2 = 0.34. The proton heating rate Q±p can be compared to 
the cascade power F that would be present at k±pp^l in "bal- 
anced" (see section |272] | AW/KAW turbulence with the same 
value of 5vp. When ci = 0.75 and C2 = 0.34, the ratio Q±p/r 
exceeds 1 /2 when Ep > Edit = 0.19, where Ep is the value of E 
for thermal protons. 

Our expression for 2±p/F (equation (i30] l) may differ from 
the value of Q±p/r in the solar wind for two main reasons. 
First, our formula for F does not take into account "imbal- 
ance" (see section 12.2b . which affects the relation between F 
and 5vp in a way that is not yet understood. Second, in true 
turbulence (as opposed to randomly phased waves), a signif- 
icant fraction of the cascade power may be dissipated in co- 
herent structures i n which the flu ctuating fields are larger than 
their rms values (iDmitruk et al.l 1200 4). Proton orbits in the 
vicinity of such structures are more stochastic than in aver- 
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Fig. 4. — The values of Ei = Svi/v^j from equation i3H\ as a function of lieliocentric distance for H+, He++, and O^'. For this figure, we assume that T_i_ = Tp 

for He^^ and and that the one-dimensional velocity power spectrum pj,' ' is « k^''^ . From bottom to top, the three curves in each plot correspond to C3 = 5/3, 
C3 = 3/2, and C3 = 1.2. 



age regions, and thus C2 may be smaller in AW/KAW turbu- 
lence than in our test-particle simulations, indicating stronger 
heating. The perpendicular heating rate is very sensitive to 
the value of C2/£p; our test-particle simulations are consistent 
with 2±p/r being exp(— C2/£p). Thus, decreasing C2 leads 
to a large increase in Q±p/r when Ep < C2- Decreasing C2 also 
decreases Eait, the value of Ep at which Q±p/r ^ 1 /2; it fol- 
lows from equations ( l24l i and ( |29] l that if Ck — 2, 5Bp/Bo — 
0.845vp/vA, and ci ~ 1, then e^rit — C2/2. 

When P -C 1, stochastic proton heating by AW/KAW tur- 
bulence at ^xPp ^ 1 increases T± much more than 7]|. In 
contrast, linear proton damping of KAWs with co <C Hp and 
^±Pp ^ 1 leads almost entirely to parallel heating, and is 
only significant when the proton thermal speed is > va; i.e., 
when Pp > 1 (O uataert.1998.) . If we assume that (nonlinear) 
stochastic heating and linear wave damping are the only dissi- 
pation mechanisms for low-frequency AW/KAW turbulence,^ 
then we arrive at the following conclusions about how the cas- 
cade power in AW/KAW turbulence is partitioned between 
parallel and perpendicular heating, and between protons and 
electrons: 

1. If Pp <C 1 and Ep ^ Ecrit, then proton heating is negligi- 
ble and electrons absorb most of the cascade power. 

2. If Pp <C 1 and gp > Eait, then parallel proton heating 
is negligible, and AW/KAW turbulence leads to a com- 
bination of electron heating and perpendicular proton 
heating. 

3- If Pp ^ 1 andEp ^Ecrit, then perpendicular proton heat- 
ing is negligible, and AW/KAW turbulence results in 
a combination of electron heating and parallel proton 
heating. 

4. If pp > 1 and Ep > Ecrit, then perpendicular proton heat- 
ing, parallel proton heating, and electron heating each 
receives an appreciable fraction of the cascade power. 

KAW turbulence at k±Pp ^ 1 fluctuates over length (time) 
scales much greater than Pe (fig"'), where Pe (Hg) is the 

^ See Markovskii and HollweJ j2002l) and lMarkovskii et alj HMS) for an 
argument against this assumption. 



thermal-electron gyroradius (cyclotron frequency). Because 
of this, an electron's magnetic moment is nearly conserved 
when it interacts with KAW turbulence at k±pp ^ 1. Elec- 
tron heating by KAW turbulence at k±Pp 1 is thus primarily 
parallel heating. On the other hand, some of the fluctuation 
energy may cascade to scales <C Pp. The way that turbulence 
is dissipated at such scales is not yet well understood. 

To determine the dependence of Ei (the value of E for ther- 
mal ions) on heliocentric distance r for different ion species in 
the fast solar wind, we adopt a simple analytic model for the 
radial profiles of the solar-wind proton density, proton temper- 
ature, and m agnetic field strength. We then apply the analyti- 
cal model of IChandran and HoUwegl (120091) . which describes 
the radial dependence of the rms amplitudes of Alfven waves 
at the outer scale Lq of the turbulence, and assume that the 

velocity power spectrum f^^' is 0= k^'^^ for Lq ^ < A;x < Pp ^ 
We find that the value of Ei for protons. Helium, and minor 
ions depends strongly on C3. However, for a fixed value of C3, 
Ei is relatively insensitive to r for 2Rq < r < I AU. 

We are not yet able to determine with precision the perpen- 
dicular heating rates of different ion species as a function of r 
because of the uncertainties in the values of C2 and C3 in the 
solar wind, and because of the large sensitivity of the heating 
rates to these quantities. However, if we assume that the value 
of C2 for protons in the solar wind is close to the value of 0.34 
in our test-particle simulations, then we arrive at the following 
two conclusions. First, perpendicular proton heating is a neg- 
ligible fraction of the turbulent cascade power in the bulk of 
the explored solar wind, in which C3 is measured to be in the 
range of 1.5 - 1.7. Second, if stochastic proton heating is im- 
portant close to the Sun, then C3 must be significantly smaller 
close to the Sun than at 1 AU. (For example, if C3 — 1.2, then 
G±p/r > 0.5 for 2R,3 <r< IOORq.) 

We find that alpha particles and minor ions undergo much 
stronger stochastic heating than protons, in large part part be- 
cause the value of Ei is larger for these ions than for protons 
at equal temperatures. Depending on the values of C2 and C3, 
the stochastic heating rate per unit volume in the solar wind 
may be larger for Helium than for protons, even though He- 
lium comprises only ^ 20% of the solar-wind mass. Figure|4] 
suggests that stochastic heating is important for alpha parti- 
cles and minor ions even if C3 is as large as 3/2, since Ei is 
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then > 0. 1 over a wide range of r. However, further investiga- 
tions into the value of eg close to the Sun and the value of C2 
for (non-random-phase) AW/KAW turbulence are needed in 
order to develop a more complete and accurate picture of 
stochastic ion heating in the solar wind. 
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APPENDIX 

LEADING-ORDER CONSERVATION OF THE FIRST ADIAB ATIC INVARIANT WHEN ^T^p ~ 1 AND E < 1 

In this appendix, we consider the interaction between ions and low-frequency, 2D = 0), electrostatic fluctuations with 
A^xP ^ 1- We assume that e ^ 1, neglect magnetic-field fluctuations, and show that the leading-order non-vanishing terms in 
dH /dt are unable to cause secular perpendicular ion heating. We set B — Bqz, where Bo is a constant. The time derivative of the 
ion's guiding-center position, defined in equation (|9|l, is then given by 

dR ^ cExz 

— =v-zH . (Al) 

dt Bo 

Since 8 <C 1, the particle's orbit in the jcy-plane during a single gyroperiod is approximately a circle of radius p = v±/Q.. We 
assume that <I> varies slowly in time, on a time scale of ^ 8^'£2^', with d<P/dz = 0. We introduce two related forms of "gyro- 
averages." First, if h is some physical property of a particle, such as its energy or guiding-center velocity, then we define the 
gyro-average of h to be 

£2 rr+Tt/a 

Ht))^^ h(ti)dh. (A2) 

271 Jt-nia 

Second, for a general function of position and time g{v,t) satisfying dg/dz = 0, we define the gyro-average of g{r,t) for particles 
with perpendicular velocity and guiding center R to be given by 

(g(r,r))R,„^ = — g{R + s{Q),t)dQ, (A3) 
where s = xp cos(9) + yp sin(9) is the vector illustrated in figure |5] 



® 
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Fig. 5. — In the small-e limit, an ion's trajectory in the xy plane is approximately a circle centered on its guiding-center position R. 



To simplify the notation, we define 

g{R,t) = {g{r,t))R.,^, (A4) 

where the functional dependence of J on is not explicitly written. If g varies slowly in time at a fixed point in space (e.g., on 
the time scale £^'£2^'), then g{R,t) is (to leading order in e) equivalent to a time average over one cyclotron period of g(r,f) 
evaluated at the position r(f ) of a particle with guiding center R: 

l(R,f) = — / g{r{ti),ti)dh. (A5) 

271 Jt-K/Q. 



13 

Thus, if we take the gyro-average of the "particle property" dR/dt in equation (lAlb using equation ( IA2l i. we find that 



dR\ ^ c 

— > = V'zZ + — 

dt / Bo 



= vzz+ — (E(r,r))R,,,^ X z. (A6) 



We consider electrostatic fluctuations with dA/dt = 0, and thus, E = — V<I>. Omitting the explicit time dependence of O and <P 
to simplify the notation, we can write the gyro-average of d<P/dx as 

/4>(r + x5)-4>(r)\ ,. 0(R + x8)-0(R) 30 
— ) = hm ( ^ ^ — ) = hm ^ ^ — = — (A7) 



3x / jj 8^0 \ 5 / R ^ 

where x is a unit vector in th e x d irection, and d/dx' denotes a partial derivative with respect to the x-component of the guiding- 
center position R. Equation (IA6b can thus be re-written as 

^WvzZ-^V'Oxz, (A8) 
dt I Bo 

where V indicates a gradient with respect to the coordinates of the guiding-center position R. Since we have assumed d/dz = 0, 

equation (lAST l impHes that {dR/dt) ■ V'O = 0. 
We now integrate equation (|8]l for an integral number of cyclotron periods, from fa to t\y = t^ +Ndt, where bt = 2n/Q. and 
> e^^. We define fo = fa + §f/2 and tj = fy_i + 5f for any integer j. Since we have assumed that dA/dt = 0, the integral of 

equation (|8]l can be written 

H{ti,)-H{t,)^qY, / -^dt. (A9) 



In analogy to equation ( IA7b . it is straightforwai'd to show that d<P/dt = {d/dt)<i>. We can thus re-write equation (|A9b as 

"3f 



t.' 30 

H{t^)-H{t,)=qY^ — (R(r,),0)5f. (AlO) 

;=o 



The time scale on which <I>(R(f),f) changes by a factor of order unity is e^'Sf. This is because <t> changes slowly in time at 
a fixed point in space, k^p ^ 1, and dR/dt ^ ev^. As a result 3<I>/3f is approximately constant within each time interval of 
duration 5f . The right-hand side of equation dAlOl ) is therefore a discrete approximation of the integral of qd<P/dt from fa to fb, 
with a fractional error of order e, so that 

H{tb)-H{t,)^q £' ^ *(R(f ) ,t)dt + ..., (All) 

where the ellipsis (. . .) represents corrections that are higher order in e. The right-hand side of equation ( lAl lb can be re-written 
in terms of the total time derivative of <I>, yielding 

rtb d — r'b dR — 

H{t^)-H{t,)=q -<P{R{t),t)dt ~ q —■V''P(R{t),t)dt + ... (A12) 

Jra dt Jt.^ dt 

Sinc e V'cS > is nearly constant during a single time interval of duration dt, the second integral on the right-hand side of equa- 
tion ( IA12b satisfies the relation 



/' 

Jts, 



^ • V'0(R(f),f)^f = E / ^dt)- V'0(R(f,),f,) + . . . (A13) 



dt j^Q \Jtj-&t/2 dt 



The integral within parentheses on the right-hand side of equation ( IA13b is equivalent to (dR/dt) 5t evaluated at f = tj. From 



equation 
equation 



(dR/dt) ■ V'O = 0. Thus, the right-han d side of equation ( IA13I ) and the second integral on the right-hand side of 
A12b vanish to leading order in £. Equation (IA12b thus becomes 

//(fb)-//(fa) -?^(R(fb),fb)-?5;(R(fa),fa) + ... (A14) 

The right-hand side of equation (IA14I) remains < ^5<t>p, regardless of how large the interval (fb — fa) becomes. Thus, to leading 
order in e, there is no secular change in the particle energy H, consistent with the near-conservation of the first adiabatic invariant 
in the small-e, small-co/il limits. 
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